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Q . Abstract 

Q I Recent results from jet production in deep inelastic ep scattering at HERA are re- 

■ viewed. The values of as{Mz) extracted from a QCD analysis of the data are presented. 

1 Introduction 

I Jet production in neutral-current (NC) deep inelastic ep scattering (DIS) provides a test of per- 
turbative QCD (pQCD) calculations and of the parametrisations of the proton parton densities 
Q ■ (PDFs). Jet cross sections allow the determination of the fundamental parameter of QCD, the 

. strong coupling constant a^, and help to constrain the parton densities in the proton. 
D ! Up to leading order (LO) in as, jet production in NC DIS proceeds via the quark-parton 
model (QPM) {Vq — > q, where V = •y or boson-gluon fusion (BGF) (Vg qq) and QCD- 
^ ' Compton (QCDC) {Vq — > qg) processes. The jet production cross section is given in pQCD by 
the convolution of the proton PDFs and the subprocess cross section. 
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(icTjet = ^ dx faix,flF) daaix,asiHR),flR,flF), 

where x is the fraction of the proton's momentum taken by the interacting parton, fa are the 
proton PDFs, /i^? is the factorisation scale, (Jq is the subprocess cross section and /i/j is the 
renormalisation scale. 

All the data accumulated from HERA and fixed-target experiments have allowed a good 
determination of the proton PDFs over a large phase space. Then, measurements of jet pro- 
duction in neutral current DIS provide accurate tests of pQCD and a determination of the 
fundamental parameter of the theory, a^. 

At high scales, calculations using the DGLAP evolution equations have been found to give a 
good description of the data up to next-to- leading order (NLO). Therefore, by fitting the data 
with these calculations, it is possible to extract accurate values of as and the gluon density of 
the proton. However, for scales of i?;^* ~ Q, where E^r^^ is the jet transverse energy and Q is the 
exchanged photon virtuality, and large values of the jet pseudorapidity, t]^'^^, large discrepancies 
between the data and the NLO calculations have been observed at low x. This could indicate 
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a breakdown of the DGLAP evolution and the onset of BFKL effects. These discrepancies can 
also be explained by assigning a partonic structure to the exchanged virtual photon or a large 
contribution of higher order effects at low . 

2 Inclusive jet cross sections 

Inclusive jet cross sections have been measured in the Breit frame using the /cy-cluster 
algorithm in the longitudinally invariant mode. The measurements were made in the kinematic 
region given by > 125 GeV^ and —0.7 < cos 7 < 0.5, where 7 is the angle of the struck 
quark in the quark-parton model in the HERA laboratory frame. The cross sections refer to 
jets of > 8 GeV and —2 < 77^* < 1.8, where i?:^^ and r^g* are the jet transverse energy 
and pseudorapidity, respectively, in the Breit frame. 

The use of inclusive jet cross sections in a QCD analysis presents several advantages: inclusive 
jet cross sections are infrared insensitive and better suited to test resummed calculations and 
the theoretical uncertainties are smaller than for dijet cross sections. 

Figure [H shows the inclusive jet cross section as a function of and -E^^g. The dots 
are the data and the error bars represent the statistical and systematic uncertainties; the 
shaded band displays the uncertainty on the absolute energy scale of the jets. The measured 
cross sections have a steep fall-off, by five (four) orders of magnitude within the measured 
(-E'tb) range. The lines are the NLO calculations using DISENT with different choices of 
the renormalisation scale {fiR = Q or -E^b)- '^^^ calculations describe reasonably well the 
and i?;^B dependence of the cross section for Q'^ > 500 GeV^ and i?;^^ > 15 GeV. At low 
and low -Ey g, the measurements are above the calculations by about 10%, which is of the same 
size as the theoretical uncertainties (see below). Therefore, for the extraction of a.^, the phase 
space was restricted to high and high -E^g. 

The experimental uncorrelated uncertainties for these cross sections are small, ~ 5%. The 
uncertainty coming from the absolute energy scale of the jets is also small, ~ 5%. The theo- 
retical uncertainties comprise 5% from the absent higher orders, 3% from the uncertainties of 
the proton PDFs and 5% from the uncertainty in the value of as{Mz) assumed. The parton- 
to-hadron corrections are 10% with an uncertainty of 1%. 

2.1 Determination of as 

The method JT] used by ZEUS to extract as exploits the dependence of the NLO calculations on 
as{Mz) through the matrix elements (a ~ A ■ ag + B ■ a"^) and the proton PDFs (^^(M^) value 
assumed in the evolution). To take into account properly this correlation, NLO calculations 
were performed using various sets of PDFs which assumed different values of as{Mz)- The 
calculations were then parametrised as a function of as{Mz) in each measured or -E^g 
region. From the measured value of the cross section as a function of in each region of Q^, 
a value of as{Mz) and its uncertainty were extracted using the parametrisations of the NLO 
calculations. 

From the inclusive jet cross section for > 500 GeV^, the value 

as{Mz) = 0.1212 ±0.0017 (stat.) t'oS (e^p.) t'omlr (th.) 

was extracted using the method explained above. The experimental uncertainties are dominated 
by the uncertainty on the absolute energy scale of the jets (1%). The theoretical uncertainties 
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Figure 1: Inclusive jet cross sections TQ as a function of (a) and (b) -E^^. 
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are: 3% from the absent higher orders, 1% from the PDFs and 0.2% from the hadronisation 
corrections. This determination is compatible with other independent extractions performed at 
HERA and with the current world average (see figure EK)- Further precision in the extraction 
of as{Mz) from inclusive jet cross sections depends upon further experimental and theoretical 
improvements. 

The QCD prediction for the energy-scale dependence of as has been tested by determining 
as from the measured differential cross sections at different scales pQ . From the measured cross 
section as a function of -E^g, in each region of -E^g, a value of ^^(ii^;^^) was extracted. The 
result, shown in figure |2b (triangles), is compatible with the running of as as predicted by 
QCD (shaded band) over a large range in the scale. Figure Eb also shows other studies of the 
energy-scale dependence of as from HERA: all the results are compatible with each other and 
with the QCD prediction. This constitutes a test of the scale dependence of as between /i = 8.4 
and 90 GeV. 



3 Parton evolution at low x 

Dijet data in DIS may be used to gain insight into the parton dynamics at low x. The evolution 
of the PDFs with the factorisation scale can be described by the DGLAP evolution equations 
which sum the leading powers of terms like as log in the region of strongly ordered transverse 
momenta /ct- This presciption describes successfully jet production at high Q^. However, the 
DGLAP approximation is expected to break down at low x since when logQ^ <C logl/x, the 
terms proportional to as log 1/x become important and need to be summed. This is done in the 
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Figure 2: (a) Summary of extracted as{Mz) values at HERA, (b) Summary of extracted 
values as a function of /i at HERA. 



BFKL evolution equations; the integration is taken over the full kx phase space of the gluons 
with no fc-r-ordering. 

Another approach, the CCFM evolution equations with angular-ordered parton emission, is 
equivalent to the BFKL approach for x — > and reproduces the DGLAP evolution equations at 
large x. Thus, the properties of the dijet system, which depend on the dynamics of the ladder, 
can be studied to determine whether the cascade has a /cT-ordered or unordered evolution. 

Deviations from the DGLAP approach at small x can be tested experimentally. At small 
X it is expected that parton emission along the exchanged gluon ladder should increase with 
decreasing x. A clear experimental signature of this effect would be that the two outgoing hard 
partons are no longer back-to-back and so an excess of events at small azimuthal separations 
should be observed. 

Values of the azimuthal separation of the two hard jets in the 7*p centre-of-mass frame, A0*, 
different than vr can occur in the DGLAP approach only when higher order contributions are 
included. On the other hand, in the BFKL and CCFM approaches, the number of events with 
A0* < vr should increase due to the partons entering the hard process with large kx- 

3.1 Azimuthal jet separation 

To test the predictions of the different approaches, dijet cross sections have been measured P] 
using the fcr-cluster algorithm in the longitudinally inclusive mode in the 7*p centre-of-mass 
frame. The measurements were made in the kinematic region given by 5 < < 100 GeV^ 
and 10"'^ < X < 10^^. The cross sections refer to jets of > 5 GeV, —1 < ?7lab < 2. 5 and 

max > GeV, where E^ is the jet transverse energy in the YP centre-of-mass frame. 

FigureEti' shows the measured dijet distribution as a function of A0*. A significant fraction of 
events is observed at a small azimuthal separation. Since a measurement of a multi-differential 
cross section as a function of x, and A0* would be very difficult due to large migrations, the 
fraction of the number of dijet events with an azimuthal separation between and a, where a 
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Figure 3: (a) Azimuthal separation between jets [21 • (b) Ratio S" as a function of Bjorken x 
and compared with predictions from NLO QCD calculations j2]. 



was taken as a = |7r, was measured instead. The fraction S, defined as 

^_ J,^N,^ot{A4>*,x,Q')dA^* 
/;iV2jet(A0*,x,Q2)rfA0*' 

is better suited to test small-x effects than a triple differential cross section. 

The measured fraction S" as a function of Bjorken x in different regions of is presented 
in figure Eb- The data rise towards low x values, especially at low Q^. The NLO predictions 
from DISENT, which contain kx effects only in the first order corrections, are several standard 
deviations below the data and show no dependence with x. On the other hand, the predictions 
of NLOJET, which contain kr effects at next-to- lowest order, provide an accurate description 
of the data at large and large x. However, they fail to describe the increase of the data 
towards low x values, especially at low Q^. 

Figure shows the data compared with the predictions of RAPGAP with direct only and 
resolved plus direct processes. A good description of the data is obtained at large and large 
X. However, there is a failure to describe the strong rise of the data towards low x, especially 
at low Q^, even when including a possible contribution from resolved virtual photon processes, 
though the description in other regions is improved. 

If the observed discrepancies are due to the influence of non-ordered parton emissions, models 
based on the color dipole or the CCFM evolution could provide a better description of the 
data. Figure shows the data compared with ARIADNE (dotted lines) and two predictions 
of CASCADE which use different sets of unintegrated parton distributions. These sets differ 
in the way the small-/cT region is treated: in Jung2003 the full splitting function, i.e. including 
the non-singular term, is used in contrast to JS2001, for which only the singular term was 
considered. The predictions of ARIADNE give a good description of the data at low x and 
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Figure 4: Ratio ^2] as a function of x and compared with predictions from RAPGAP (a) 
and ARIADNE and CASCADE (b). 



Q^, but fail to describe the data at high Q'^. The predictions of CASCADE using JS2001 he 
significantly above the data in all x and regions, whereas those using Jung2003 are closer 
to the data. Therefore, the measurement of the fraction S is sensitive to the details of the 
unintegrated parton distributions. 



4 Internal structure of jets 

The investigation of the internal structure of jets gives insight into the transition between 
a parton produced in a hard process and the experimentally observable jet of hadrons. The 
internal structure of a jet depends mainly on the type of primary parton from which it originated 
and to a lesser extent on the particular hard scattering process. QCD predicts that at sufficiently 
high -E;^*, where fragmentation effects become negligible, the jet structure is driven by gluon 
emission off the primary parton and is then calculable in pQCD. The lowest non-trivial order 
contribution to the jet substructure is given by order aag calculations. 

The internal structure of the jets can be studied by means of the mean subjet multiplicity. 
Subjets are resolved within a jet by reapplying the algorithm on all particles belonging to 
the jet until for every pair of particles the quantity dij = mm{ET,i, EtjY ■ {{rji—rjjY + [ipi — ipjY) 
is above dcut = Umt ■ {Ej^Y. All remaining clusters are called subjets. The subjet structure 
depends upon the value chosen for the resolution parameter Ucnt- 

The mean subjet multiplicity has been measured [2] for jets using the algorithm in 
the HERA laboratory frame with i?;^* above 15 GeV and — 1 < rf^^ < 2, in the kinematic 
range given by > 125 GeV^. Figure Et- shows the mean subjet multiphcity for a fixed 
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value of i/cnt of 10~^ as a function of i?;^*. It decreases as i?;^* increases, i.e. the jets become 
more collimated. The experimental uncertainties are small (fragmentation model uncertainty 
< 1%, the uncertainty on the absolute energy scale of the jets is negligible). The detector and 
hadronisation corrections are < 10% and < 17%, respectively, for > 25 GeV. The data 
are compared to the LO and NLO predictions of DISENT. The LO calculation fails to describe 
the data, whereas the NLO calculations provide a good description. These measurements are 
sensitive to and have been used to extract a value of as{Mz). The result is 

as{Mz) = 0.1187 ±0.0017 (stat.) t'oZt (exp.) Toow (th-)- 
This value is compatible with the world average and with previous measurements (see figure|2ti')- 
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Figure 5: (a) Mean subjet multiplicity fP as a function of i?;^*. (b) 1 — cr contours in the (a^, o.^) 
plane 



5 Event shapes 

A complementary extraction of a^, also using the details of the hadronic final state in DIS, 
comes from the study of the event shape variables, like thrust or jet broadening. Event shape 
variables are particularly sensitive to the details of the non-perturbative effects of hadronisation 
and can be used to test the models for these effects. Recently, new developments with regard to 
power-law corrections have prompted a revived interest in the understanding of hadronisation 
from first principles. In this type of analysis, the data are compared to model predictions which 
combine NLO calculations and the theoretical expectations of the power corrections, which are 
characterised by an effective coupling ao- Previous results supported the concept of power 
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corrections in the approach of Dokshitzer et al. but a large spread of the results suggested that 
higher order corrections were needed. Now, resummed NLL calculations matched to NLO are 
available and so it is possible to study event shape distributions instead of only their mean 
values. 

Event shape distributions (thrust, jet broadening, the jet mass p and the C parameter) 
have been measured |3] for particles in the current hemisphere in the kinematic region given 
by 14 < Q < 200 GeV and 0.1 < y < 0.7, where y is the inelasticity variable. Predictions 
consisting of NLO calculations using DISASTER++, resummed calculations matched to NLO 
and power corrections have been fitted to the data, leaving and as free parameters. A 
good description of the data by the predictions was obtained at high Q^, though the description 
at low was poorer. Figure shows the la-contour results from the fit in the Us-otQ plane. 
The spread observed in previous studies is much reduced when the resummed calculations are 
included. A clear anti-correlation between «<, and is found for all variables. A universal 
value for ao of 0.5 at the 10% level was obtained, in agreement with the previous results, but 
with a smaller spread. There is still a sizeable theoretical uncertainty for both as and ao, of 
the order of 5%, which is as large as the experimental uncertainties. This uncertainty comes 
from the absent higher order corrections. 

6 Conclusions 

HERA has become a unique QCD-testing machine due to the fact that at large scales consid- 
erable progress in understanding and reducing the experimental and theoretical uncertainties 
has led to very precise measurements of the fundamental parameter of the theory, the strong 
coupling constant a^. The use of observables resulting from jet algorithms leads now to deter- 
minations that are as precise as those coming from more inclusive measurements, such as from 
r decays. To obtain even better accuracy in the determination of QCD, further improvements 
in the QCD calculations are needed, e.g. next-to-next-to-leading-order corrections. 

At low values of x and Q^, considerable progress has also been obtained in understanding the 
mechanisms of parton emission, though the interplay between the DGLAP, BFKL and CCFM 
evolution schemes has still to be fully worked out. Further progress in this respect needs both 
more experimental and more theoretical work. 
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